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Abstract: There a sizable and systematic discrepancy between experimental data on 
the bb production in pp, 7p and 77 collisions and existing theoretical calculations within 
perturbative QCD. Before interpreting this discrepancy as a signal of new physics, it is im- 
portant to understand quantitatively the ambiguities of conventional calculations. In this 
paper the uncertainty coming from renormalization and factorization scale dependence 
of finite order perturbation calculations of the total cross section of bb production in pp 
collisions is discussed in detail. It is shown that the mentioned discrepancy is reduced 
significantly if these scales are fixed via the Principle of Minimal Sensitivity. 
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1. Introduction 



Heavy quark production in hard collisions of hadrons, leptons and photons has been 
considered as a clean test of perturbative QCD. It has therefore come as a surprise that 
the data on bb production in pp collisions at the Tevatron [1, 2], 7p collisions at HERA 
[3, 4] and 77 collisions at LEP2 [5, 6], some of them reproduced in Fig. 1, lie systematically 
by a factor of about 2-4 above the median of current theoretical calculations. For the first 
process the data are sufficiently copious to measure also the differential distribution in 
transverse momentum of the produced b quark. Although the data have sizable errors, this 
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Figure 1: Data on bb production in pp [1], 7p [3, 4] and 77 [5, 6] collisions. 



excess is too big to be accommodated within the current calculations. The comprehensive 
review of current status and related problems of the comparison of data on heavy quark 
production in all three collisions with QCD calculations can be found in [7, 8, 9]. 

The discrepancy between the DO and CDF data on bb production and theoretical cal- 
culations [10], shown in Fig. la, has led to suggestions that it might represent the mani- 
festation of effects of unintegrated gluon distribution function within the rCT-factorization 
approach [11, 12], or even a signal of supersymmetry [13, 14]. The agreement between 
the Tevatron data and the calculations [11, 14], reproduced in Fig. 2, is, indeed, quite 
impressive. On the other hand, before invoking new physics it is important to check that 
the data cannot be explained by more sophisticated application of QCD phenomenology. 

The QCD calculations of heavy quark production in hard collisions of hadrons, leptons 
and photons depend on a number of inputs: a s , parton distribution functions (PDF) 
of colliding hadrons or photons, fragmentation functions of b quarks, masses of heavy 
quarks, and the choice of renormalization (RS) and factorization (FS) scales /x and M. 
Quite recently, proper parameterization of fragmentation function has been argued [9] to 
be one of the possible sources of the disagreement of the CDF data [2] with NLO QCD 
calculations. There are also attempts to go beyond fixed order perturbative calculations 
by resumming the effects of large logs of the type ln(S'/4m^) [15] or \n(p T /m b ) [16]. 
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Figure 2: DO and CDF data on bb production compared to calculations using the k~T- 
factorization approach as implemented in CASCADE event generator [11] (left) and the su- 
persymmetry [?] (right). 



In this paper I investigate in detail the dependence of existing fixed order (LO and 
NLO) QCD calculations [17, 18] of the total cross section <J tot (bb) in pp collisions on the 
choices of the renormalization and factorization scales. Similar analysis of bb production 
in 7p and 77 collisions will be presented elsewhere 1 . 

The paper is organized as follows. In the next Section basic facts relevant for further 
discussion are collected. In Section 3 some of the methods for selecting the renormaliza- 
tion and factorization scales are recalled and their relevance for the process under study 
discussed. This is followed in Section 4 by the discussion of the general form of the renor- 
malization and factorization scale dependence of Utotibb) at the NLO. Numerical results 
relevant for energies up to the Tevatron range are presented in Section 5 and conclusions 
drawn in Section 6. 



2. Basic facts and formulae 

The basic quantity of perturbative QCD calculations, the renormalized color coupling 
ce s (n), depends on the renormalization scale /1 in a way governed by the equation 

where for rif massless quarks (3$ = 11 — 2n//3 and f3\ = 102 — 38n//3. The solutions 
of (2.1) depend beside fj, also on the renormalization scheme (RS). At the NLO (i.e. 
taking into account first two terms in (2.1)) this RS can be specified, for instance, via the 
parameter Ars, corresponding to the renormalization scale for which a s diverges 2 . The 



1 For discussion of the specific features of bb production in 77 collisions see [19] 

2 At higher orders the /3-function coefficients f3k,k>2 may be used in addition to A^g to fix the RS. 
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coupling is then given as the solution of the equation 

4tt V a rs/ a.M l + ca s (/i)' 

where c = /3i/(47r/3 )- At the NLO the coupling a s is thus a function of the ratio /i/A RS 
and the variation of the RS for fixed scale ji is therefore equivalent to the variation of 
fi for fixed RS. To vary both the renormalization scale and scheme is legitimate, but 
redundant. Let me emphasize that the choice of the RS is as important as the choice 
of renormalization scale. The fact that the coupling a s (fi), as well the coefficients of 
perturbation expansions, depend for fixed \x on the RS also implies that the existence of 
a "natural" physical scale in the problem, like the masses of heavy quarks in our case, is 
actually of no help in specifying unambiguously the NLO calculation. I will come back 
to this point in the next Section. If not stated otherwise, I will work in the conventional 
MS RS and vary the renormalization scale \i only. 

For hadrons the factorization scale dependence of PDF is determined by the system 
of evolution equations for quark singlet, nonsinglet and gluon distribution functions 



dS(M) 

dlnM 2 
dG(M) 

din M 2 

dq m (M) 

dlnM 2 



P qq {M) ® E(M) + P qG {M) ® G(M), (2.3) 
P Gq {M)®Y,{M) + P GG {M)®G{M), (2.4) 
P NS (M) ®g NS (M), (2.5) 



where 



E(x, M) = J2 (fcfo M ) + Ux, M)) , (2.6) 
i=\ 

q NS!i (x, M) = (q t (x, M) + q^x, M)) - — S(x, M), Vi (2.7) 

n f 

The splitting functions admit expansion in powers of a s (M) 

Plj(x ,M) = W+ (f^)) 2 ^> W + .., (2.8) 

where P^\x) are unique, whereas all higher order splitting functions Pjf[\j > 1 depend 
on the choice of the factorization scheme (FS). Conversely, they can be taken as defining 
the FS, similarly as the higher order /5-function coefficients [3i,i > 2 in (2.1) define, 
together with A^g, the renormalization scheme. The equations (2.3-2.5) can be recast 
into evolution equations for qt(x, M), q^x, M) and G(x, M). 

In all calculations of this exploratory study I took rif — 4 and solved the RG equa- 
tion (2.2) numerically in order to guarantee correct behaviour of its solution for small 
renormalization scales. 
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3. Renormalization and factorization scales and their choice 



Let us first recall the situation for processes involving the renormalization scale only. For 
them the QCD contribution to corresponding cross sections has the form (a s = a s /n) 

r{Q) = a k s (p, RS) (r (Q) + n(Q/fi, RS)a s (p, RS) + r 2 (Q//i, RS)a^, RS) • • •) . (3.1) 

As an example, take the familiar ratio 

r( r- r u T + ^ Irons) = 3 

T(r -> UrfjL V^) 

where the QCD correction r T is given by the expression (3.1) with k = 1 and r = 1. 
Theoretical consistency of the expansion (3.1) implies 

JL c\ r (N) 

r^^r^ai => ^- = 0(a^), (3.3) 

j=k ^ 

i.e. the derivative of the N-th order approximant r* 7 ^ with respect to In /i is of higher 
order than the appoximant itself. This requirement determines the dependence of 
the coefficients r\ on For r\ we get, for instance, 

r^Q/pi, RS) = /cfeln ^ + n(l, RS) = kbla p(Q), (3.4) 

V A RS 

where p{Q) = /c6 ln(Q/A RS ) — rx(l,RS) is a renormalization scale and scheme invariant. 

There are two general approaches to selecting the renormalization scale applicable 
in such circumstances: the Principle of Minimal Sensitivity [20] (PMS) and the method 
based on the concept of Effective Charges [21] (EC). In the first approach the emphasize 
is put on the stability of the results with respect to variations of [i. In the absence of 
information on higher order perturbative terms in (3.1) the PMS approach is natural as 
it selects the point \x where the truncated perturbation expansion is most stable and has 
thus locally the property possessed by the all order result globally. The EC approach 
is based on the criterion of apparent convergence of perturbation expansions and the 
renormalization scale \i (and at higher orders other free parameters) is chosen in such 
a way that all higher order contributions vanish, i.e. demanding = for all N. 
Whatever the choice of the renormalization scale fi scale we make, it is, however, certainly 
useful to investigate the stability of the results with respect to the variations of fi. 

In a fixed RS, both PMS and EC approaches can be phrased in geometrical terms 
as methods of selecting on a curve defined by the expression r^(/i) some preferred 
point. The PMS approach selects the stationary point(s), whereas the EC method selects 
the point(s) given by the intersection of A N \fi) with the curve corresponding to the 
lowest nontrivial order. In Fig. 3 the \i dependence of the first three approximants 
ri N \ N = 1, 2, 3 for the specific case of the quantity r T defined in (3.2) is plotted in MS 
RS taking = 0.31 GeV. Several features of this figure are worth noting. 
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• The conventional NLO approximation in the MS RS, corresponding to the choice 
^ms = m r — 1-8 GeV, gives substantially smaller result than that corresponding to 
the stationary point. 

• Working in other RS, the same choice /x RS = m T = 1.8 GeV would give a different 
result for r T (m T ). In fact, any point on the NLO and NNLO curves can be obtained 
by setting /ij^g = m T in appropriately chosen RS! The usual argument for using the 
MS RS is that the associated a s incorporates ln47r — 7^, the artefact of dimensional 
regularization, and thus leads to smaller expansion coefficients than those of MS 
RS. But if the rate of convergence of perturbation expansions of physical quantities 
is the main criterion for the appropriate choice of the RS, than that based on the 
EC approach, rather than the MS should be selected. 

• The NLO and NNLO approximants and exhibit maxima, defining the PMS 
scales /i^ S and the corresponding PMS optimized results r^} s (Q). No physical 
meaning should, however, be attached to numerical values of fJ^ s as in different 
RS the same maximum occurs at different scale 3 . 

• At both NLO and NNLO the PMS scales fi^ s turn out to be quite close the scales 

chosen by the EC approach. This fact holds, with some modification, for more 
complicated quantities as well. 

• Although the preferred scales fJ^ s and /J^ depend on the chosen RS, the PMS and 
EC optimized results r^^Q) and r^(Q) do not! This fact follows directly from 
the geometrical meaning of the preferred points and is crucial for the uniqueness of 
the obtained predictions. 

3 For instance, M ™S = 0.38/z™ s , whereas = 2.17m™ s . 
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At higher orders and for more complicated physical quantities the situation is, however, 
not always so simple as sketched above. In some cases there is no genuine stationary point 
in the physical region and one must then look for some other measure of "most stable" 
region, or, on the contrary, there may be more stationary points to choose from. For 
physical quantities with two scales the general criterion of the EC approach does not lead 
to a unique point, but rather the curve /i(M) along which the higher order contributions 
vanish. However, even in these circumstances the PMS and EC approaches provide useful 
guidance to fixing the scales and thus defining the QCD predictions. 

I now come to the central aspect of the choice of the renormalization and factorization 
scales. As emphasized by Politzer [22] , there is no compelling reason for identifying these 
two scales. The point is that whereas the renormalization scale \x defines, roughly speak- 
ing, the lower limit on the virtualities of loop particles included in the definition of the 
renormalized coupling, the factorization scale M specifies the upper limit on virtualities 
of partons included in the definition of dressed PDF. In other words, the renormalization 
scale reflects ambiguity in the treatment of short distances, whereas the factorization scale 
comes from similar ambiguity concerning large distances. It is therefore natural to keep 
the renormalization and factorization scales as independent free parameters of any finite 
order perturbative approximations. The calculations presented in the next two Sections 
demonstrate that at the NLO the cross section a tot (bb; M , /i) depends, indeed, on these 
two scales in quite different way. 

Identifying these two scales may actually lead to misleading conclusions concerning 
the stability of perturbative expressions. It may act in both ways, faking the "stability" 
where there is is fact none, but also hiding the genuine stability region when it lies away 
from the diagonal point /j, — M. At the NLO there is an additional degree of freedom 
(in fact an infinite number of them) related to the choice of the factorization scheme, 
specified by the higher order splitting functions ■ I will not exploit this freedom and 
throughout this paper stay within the conventional MS FS. 

4. General form of a to t(QQ) 

According to the factorization theorem the total cross section of heavy quark pair pro- 
duction in pp collisions at the center of mass energy \fS has the form 

°tot{w - QQ, S) = JJdxdy Df(x, M)D?(y, M)a t ,(s = xyS, M), (4.1) 

where partonic cross sections ^(s, M) = a(ij — > bb, s, M), as well as PDF of the beam 
particles, D?(M), D P (M), depend on the factorization scale M. The theoretical calcula- 
tions of heavy quark productions in refs. [7, 17, 18] are given in terms of the heavy quark 
pole mass, i.e. uiq appearing in the expressions for is just a fixed number. 

The expression (4.1) is basically of nonperturbative nature. Fixed order perturbation 
theory enters if we insert into (4.1) the solutions of the evolution equations (2.3- 2.5) 
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with the splitting functions Pij truncated to finite order (2.8) and calculate partonic cross 
sections (Jij(s, M) as power expansion in the coupling a s (fj) 

a^s, M) = a?MaW(s) + a 3 s (^f(s, M, //) + ••• , (4.2) 

taken at the renormalization scale //, which, as emphasized above, is in general unrelated 
to the factorization scale M. Summed to all orders of ot 8 , the r.h.s of (4.2) is actually 
independent of //, whereas even the all-order sum of (4.2) depends on the factorization 
scale M. The latter dependence is cancelled by that of the PDF provided the splitting 
functions P^ in the evolution equations (2.3-2.5) are taken to all orders. This difference 
reflects the different roles played by the renormalization and factorization scales in ex- 
pressions for physical quantities. If only finite order approximantions in (4.2) and (2.8) 
are used, as is in practice always the case, the partonic cross sections aij(s,M) become 
also yU-dependent and the expression (4.1) thus becomes a function of both \x and M. 

At the NLO, i.e. taking into account the first two terms in (4.2) and (2.8), we get 
(dropping for brevity the dependence on S and specification of bb final state) 



^°(M,//) = a 2 s (fi) | JJ dxdyY^ qi {x,M) qi {y,M) (xy) + a s (^ (xy, M, //) 

2 jj dxdyZix, M)G(y, M)a s (^a^(xy, M)+ (4.3) 
jj dxdyG(x, M)G(y, M) [a^ G (xy) + a s (^ G (xy, M, y) 



where the sum is understood to run over rif quarks and antiquarks, and the relation 
between PDF of protons and antiprotons was taken into account. 

Theoretical consistency of (4.1) implies, similarly to (3.3), that the derivative of any 
finite order approximation thereof with respect to In M 2 is of higher order in a s than such 
approximation itself. Performing this derivative for (4.3) we get 

da^°(M,/i) jj dxdyG ^ M)GiyjM)WGGixy ^ M ^ )+ (44) 

2n f 

qi (x, M) qi {y, M)W qq (xy, M, //) + S(x, M)G(y, M)W qG (xy, M, //) 



dlnM 2 

r2?; , 

jj dxdy 



.i=i 

Using (2.3-2.5) and denoting / = d//dlnM 2 , the functions are given as 

+jdzPW(z)a%Uxz)} + ... (4.5) 
W qq (x, M,„) = ^{ 2k&<£> (x)+2j izPW {z)a% {xz ) J + . . . (4.6) 
W qG (x, M, fj) = {^(s) + Jdz[pV>(z)v$(xz) + P%>(z)v™(xzj\}+ ■ -(4.7) 
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Figure 4: Graphical representation of cr\^ ( M, /i) for y/S = 630 GeV: a) a two dimensional 
surface, together with its both projections, plotted for several values of fi in b) and Mine). 

Theoretical consistency of (4.3) requires that the expressions standing in (4.5-4.7) by 
vanish which, indeed, they do. For practical calculations it is suitable to rewrite the 
expression (4.1) as an integral 

VtotQQQ, S) = [ dz ®iA z , M)a i:j (s = zS, M) (4.8) 

a 

over the product of parton fluxes 

M) = J J dxdyDf(x, M)D p j (y, M)S(z - xy) (4.9) 

and partonic cross sections a{j(z,M), both of which depend on the product of fractions 
z = xy only. Because the expressions for a\j [s, M) as given in [17] correspond to \i = M, 
I have restored their separate dependence on /x and M by adding to cr^(s, M) the term 
2/5 ag ) ln(/i 2 /M 2 ). 

5. Results 

In Fig. 4 the LO approximation a^(M, /i) i.e. the sum of the terms proportional to o? s in 
(4.3), is plotted for = 630 GeV two-dimensional surface a^(M,fi), as well cLS cL 
function of M for fixed /i, and as a function of \i for fixed M. There is no stability region 
of 

<T io?(^/ i ) 111 t ne sense of a stationary point on a two-dimensional surface. Keeping 
one of the scales fixed and varying the other (solid curves in Figs. 4b-c) leads in both 
cases to monotonous dependence on the latter. If, however, we set /i = M, which means 
cutting the surface in Fig. 4a along the diagonal line, we obtain the dashed curve shown 
in both projections in the Figs. 4b-c, which has a spurious stationary point at about 
/j, = M = 1.5 GeV, illustrating the point made at the end of Section 3. 
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The NLO approximant af^°{M, //) is represented graphically in Fig. 5 for \fS = 62 
GeV, one of the energies investigated in [18], m h = 5 GeV and using GRV HO PDF. 
Setting /i = M we get the curves in Fig. 5a, where the individual contributions of 
gluon-gluon, quark- ant iquark and quark-gluon channels are plotted as well. The quark- 
antiquark channel dominates for M > 2.5 GeV, but the gluon-gluon one contributes 
significantly, too. At small /i = M the quark-quark contribution turns negative, as does 
eventually the gluon-gluon one, whereas that of the quark-gluon channel rises. The latter 
channel gives negative contribution at large values of M. To quantify the effects of the 
conventional assumption /i = M, afJ^°(M, ^ = k,M) is plotted in Fig. 5b for three values 
of k = 0.5, 1, 2. The position of the stationary point as well as the corresponding value 
of <r^°(M, /i = kM), depend on the value of n, though not dramatically. Note that 
the positions of the maxima in Fig. 5b occur close to the intersections of the LO and 
NLO curves, defining the EC predictions. The spread between the three NLO curves at 
M = nib may be used as one of the measures of theoretical uncertainty of the conventional 
NLO calculations. The standard method of estimating this uncertainty is to set M = fi 
and vary this common scale within the interval (mfe/2,2mb). Note, however, that the 
resulting band of values of af^ (Ki7ib, nrrib) depends on the choice of RS! 
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Keeping \x and M as independent parameters and representing crf^ (M, /x) as two- 
dimensional surface in Fig. 5c or by the corresponding contour plot in Fig. 5d, reveals 
the saddle point located approximately at il — 2.4 GeV, M = 3 GeV. This saddle point, 
which defines the PMS result, lies close to the curve in Fig. 5d along which ofjf (M, /i) = 
af^(M, //), defining the results preferred by the EC approach. Contrary to the PMS the 
Effective Charges approach does not lead to a unique prediction, but the whole range of 
values. At y/S = 62 GeV the saddle point lies also close to the point fx = M = m b /2 = 2.5 
GeV, which is conventionally taken as the lower limit on the variation of the common scale 
fj, — M — m b . As the collision energy increases, the situation changes in two respects: 




M [GeV] M 



Figure 6: The same as in Fig. 5a,d but for \/~S = 630 GeV. 

• The relative contribution of the quark-antiquark channel decreases, its dominance at 
low energies being taken over by that of the gluon-gluon channel. The gluon-quark 
channel grows in importance at low values of M, roughly M < 2 GeV. 

• The position of the saddle point moves in a highly "nondiagonal" manner away from 
the line il — M and thus away from the conventional choice il = M. However, at 
all energies the saddle point stays close to the curve defining the EC results. 

These changes are illustrated in Figs. 6, which shows the same plots as those in Fig. 
5a,d but for = 630 GeV. The energy dependence of the position of the saddle point 
(Hsad(S) : M sad (S)), represented by solid curves in Fig. 7a, illustrates yet again the different 
role of renormalization and factorization scales. Whereas M sa( i{S) stays constant up to 
about \/S ~ 200 GeV and then rises roughly as some power of S, n sa d(S) first dips down 
to Usad ~ 0.8 GeV at about v^S* ~ 500 GeV, before starting to rise at about the same 
rate as M sac i(S) but remaining far below it. This large disparity between the values of 
renormalization and factorization scales at the saddle point provides the main reason for 
keeping these two scale as independent free parameters. 

The "nondiagonal" position of the saddle point at higher energies leads to increasing 
difference between the PMS-optimized and conventional MS results for af ^°(S). This 
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is illustrated in Fig. 7b which shows that the former rises more steeply in the region 
100 ;$ \A ~ 500 GeV than the standard calculations, after which they start to approach 
them again. To quantify this difference I plot in Fig. 7c the ratio: 



R PMS (S) 



a 



NLO 

tot 



(M sac i, Li S ad) 



a 



NLO 

tot 



(m b ,m b , MS) 



together with the ratia 



R K (S) 



a 



NLO, 
tot 1 



KWlj, Km b , MS) 



a 



NLO 

tot 



(m b ,m b ,MS) 



(5.1) 



(5.2) 



for k = 0.5 and k = 2. The PMS-optimized scales thus lead to predictions that are 
higher than those evaluated with the conventional choice fj, = M = m b . For low energies 
the former are close to those corresponding to conventional calculations with k = 0.5, 
but at high energies R PMS (S) rises steeply, reaching its maximum R PMS (S) = 2.2 at 
about = 400 GeV. At ^fS = 1800 GeV R PMS (S) ~ 1.9, which is about the same 
number as the factor characterizing the excess of DO data over the conventional NLO QCD 
predictions in MS, which assume fj, = M = m b . The fact that the difference between the 
results based on PMS-optimized and conventional choices of scales is a sensitive function 
of the collision energy provides a way for checking the phenomenological relevance of the 
optimization procedure based on the PMS idea of local stability. 
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Figure 7: a) Energy dependence of the position of the saddle point of (bb) for m b = 5 
GeV using GRV (solid curves) and MRS (dashed curves) PDF of the proton. The dotted curves 
show the same but for m b = 11 GeV. The dash-dotted curve shows the energy dependence of 
the average pt of the 6-quark. b) The energy dependence of (bb) in the conventional MS 
RS (dotted curves) as well as the PMS-optimized prediction (solid curves), for two values of 
m b . c) The ratia i? 1 / 2 , R 2 and R PMS for m b = 5 GeV using the GRV (solid curves) and MRS 
(dashed curve) PDF. The dotted curves correspond to GRV PDF and m b = 11 GeV. 



The results presented above have been obtained using fixed order perturbative QCD 
and could therefore be modified by the inclusion of the resummation of small-x and/or 
large pt logs. Note that at the Tevatron the momentum fractions carried by partons 
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producing the bb pair are typically larger than 0.006, which is where the region of "small- 
x" usually starts. So one might expect the approach based on unintegrated PDF as used 
in CASCADE to improve the agreement with data, though it is not obvious by how much 
as CASCADE uses only a LO hard scattering cross sections. On the other hand, R PMS 
peaks at — 400 GeV, for which the mean parton fractions (x) > 0.025 are safely 
outside the small-a; region. Consequently, in the energy range where it is largest the 
difference between standard NLO calculations in MS RS and PMS-optimized ones will 
likely survive the effects of small-x resummation. 

The LO calculation of the energy dependence of mean p? of the produced 6-quarks, 
displayed by the dash-dotted curve in Fig. 7a, shows that it rises only very slowly and 
approaches 7 GeV at the Tevatron range. This indicates that the resummation of large 
logs ln(pT/mb) may be important for the tails of pr distributions, but cannot be expected 
to influence significantly the results for the total cross section. 

The results discussed so far were obtained with the GRV HO parameterization of PDF 
of the proton. To check the dependence of the above conclusions on the choice of PDF, 
all the calculations were redone also with the MRS(R2) PDF used in [1]. As illustrated 
in Fig. 7a, there is noticeable difference in the location of saddle point at low energy, 
where the GRV saddle points are closer to the "diagonal" ones than the MRS(R2) ones, 
but as the energy increase both parameterizations give essentially the same results not 
only for the for the positions of the saddle points, but also for the energy dependence of 
the corresponding optimized total cross sections. This is illustrated in Fig. 7c, where the 
ratio R PMS is plotted for both parameterizations. The conventional ratia i? 1//2 ,i? 2 are so 
close for both parameterizations, that only those corresponding to GRV are displayed. 

To get some feeling for the dependence of the above conclusions on the mass of the 
6-quark, Figure 7 shows also the results obtained for m b = 11 GeV, which is roughly 
the transverse mass ttit = \J ml + p\ at the Tevatron. Clearly, the difference between 
standard and PMS-optimized calculations decreases with increasing mj and the region of 
largest R PMS shifts, as expected, to larger energies. 

6. Summary and conclusions 

The results of this exploratory study demonstrate the sensitivity of the NLO QCD cal- 
culations of the total cross section af ^°(bb) to the choice of the renormalization and 
factorization scales. In particular they show that the conventional assumption of identi- 
fying these two scales leads at high energies far away from the region of local stability as 
quantified by the Principle of Minimal Sensitivity criterion. As a result, the conventional 
NLO calculations yield results that are lower than the PMS-optimized ones, at the Teva- 
tron energies by a factor of about 2, which is approximately the factor needed to bring 
the conventional calculations into agreement with the DO and CDF data. 

This observation cannot, however, be directly used as an explanation of the excess 
of the Tevatron data over conventional NLO predictions because it concerns the total 
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cross section rather than the differential distribution da/dpr measured by the DO and 
CDF Collaborations. The ongoing analysis [23] of renormalization and factorization scale 
dependence of such differential distribution indicates that the situation for such distribu- 
tion is similar to that discussed in this paper, with additional dependence of the stability 
region on pt of the 6-quark. This paper will also address the question of the stability of 
the NLO QCD calculations of the top quark production at both Tevatron and LHC. 

I also don't claim that the choice of the renormalization and factorization scale via the 
PMS-optimization explains fully the observed discrepancy between data and conventional 
NLO QCD calculation of afjf (bb) or that it is the only effect that might be responsible 
for such discrepancy. The purpose of this paper is rather to point out yet another source of 
theoretical uncertainty of higher order QCD calculations of quantities like afj^°(bb) that 
should be properly taken into account before drawing conclusions from the comparisons 
with experimental data. 

Acknowledgments 

This works has been supported by the Ministry of Education of the Czech Republic under 
the project LN00A006. 

References 

[1] B. Abbott et al. (DO Collab.): Phys. Rev. Lett. 85 (2000), 5068, (hep-ex/0008021) 

[2] D. Acosta et al. (CDF Collab.): Phys. Rev. D65 (2002), 052005, (hep-ex/0111359) 

[3] C. Adloff et al (HI Collab.): Phys. Lett. B467 (1999), 156 and Erratum: Phys. Lett. B518 
(2001), 331 

[4] J. Breitweg et al. (ZEUS Collab.): Eur. Phys. J. C18 (2001), 625, (hep-ex/0011081) 
[5] M. Acciarri et al. (L3 Collab.): Phys. Lett. B503 (2001), 10 
[6] OPAL Collab.: OPAL Physics Note PN455, August 2000 

[7] S. Frixione: Proceedings of International Europhysics Conference on HEP, Budapest, July 
2001, hep-ph/0111368 

[8] P. Nason: Proceedings of Interbational Symposium on Lepton and Photon Interactions at 
High Energy, Rome, July 2001, hep-ph/01 11024 

[9] M. Cacciari and P. Nason: Phys. Rev. Lett. 89 (2002), 122003-1 

[10] S. Frixione and M. Mangano: Nucl. Phys. B483 (1997), 321, (hep-ph/9605270) 

[11] H. Jung: Phys. Rev. D65 (2002), 034015 

[12] A. V. Lipatov, V. A. Saleev and N. P. Zotov: hep-ph/01 121 14 



13 



[13] E. Berger et al: Phys. Rev. Lett. 86 (2001), 4231, (hep-ph/0012001) 
[14] E. Berger: hep-ph/0201229 

[15] S. Catani, M. Ciafaloni and F. Hautman: Nucl. Phys. B366 (1991), 135 

[16] M. Cacciari, M. Greco and P. Nason: JHEP 05 (1998), 007, (hep-ph/9803400) 

[17] P. Nason, S. Dawson and R. K. Ellis: Nucl. Phys. B303 (1988), 607 

[18] G. Altarelli, M. Diemoz, G. Martinelli and P. Nason: Nucl. Phys. B308 (1988), 724 

[19] J. Chyla, hep-ph/0111469. 

[20] P. M. Stevenson, Phys. Rev. D21 (1981), 2916 

[21] G. Grunberg, Phys. Rev. D29 (1984), 2315 

[22] H. D. Politzer, Nucl. Phys. B192 (1984), 493 

[23] J. Srbek: in preparation 



14 



